Adsorption of para-Hydrogen on Pullerenes 
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Adsorption of para-hydrogen on the outer surface of a single fullerene is studied theoretically, by 
means of ground state Quantum Monte Carlo simulations. We compute energetics and radial density 
profiles of para-hydrogen for various coverages on a variety of small fullerenes. The equilibrium 
adsorbed solid monolayer is commensurate with the surface of the fullerene; as the chemical potential 
is increased, a discontinuous change is generally observed, to an incommensurate, compressible layer. 
Quantum exchanges of hydrogen molecules are absent in these systems. 

PACS numbers: 67., 68.08.De, 68.43.-h, 68.60.-p, 68.65.-k, 81.05 Tp 



I. INTRODUCTION 

Low temperature adsorption of highly quantal fluids, 
such as helium or para- hydrogen (p-H 2 ) on the outer sur- 
face of a fullerene ( "buckyball" ) can provide insight into 
physical properties of a quantum many-body system con- 
fined to spatial regions of nanometer size. As the diam- 
eter of the fullerene is increased, the properties of the 
adsorbatc ought to interpolate between those of a cluster 
with a solvated impurity, and those of an adsorbed film 
on an infinite substrate. 

In this paper, we consider adsorption of p-H 2 on a sin- 
gle fullerene C;, with 1=20, 36, 60 and 80. All of these 
molecules are strong adsorbers, and very nearly spherical. 
Background for this study is provided by the wealth of 
theoretica l 1 ! 2 !?! 4 '^ and experimenta l?!?!?! 1 '' 1 ' 11 ! 12 ! 1 ? work, 
spanning over two decades, aimed at investigating the 
properties of adsorbed P-H2 films on various substrates. 
This work is also inspired by recent theoretical results on 
adsorption of helium on buckyballs^ 4 * 1 ^ 

A fluid of J9-H2 molecules is an interesting physical sys- 
tem for a number of reasons. Because a P-H2 molecule is 
half as light as a helium atom, zero-point motion can be 
expected to be quite significant; each molecule is a spin- 
zero boson, and therefore it is conceivable that, at low 
enough temperature, a P-H2 fluid might display phys- 
ical behavior similar to that of fluid helium, including 
superfluidity^ Unlike helium, though, bulk P-H2 solidi- 
fies at low temperature (T c w 14 K); this prevents the 
observation of phenomena such as Bosc Condensation 
and, possibly, superfluidity, which are speculated to oc- 
cur in the liquid phase below T w 6 K. Solidification is 
due to the depth of the attractive well of the potential 
between two hydrogen molecules, which is significantly 
greater than that between two helium atoms. Several, 
attempts have been madeiii&iiLSSi to supercool bulk liq- 
uid P-H2, but the search for superfluidity (in the bulk) 
has so far not met with success. 

Confinement, and reduction of dimensionality, are 
widely regarded as plausible avenues to the stabilization 
of a liquid phase of p-^2 at temperatures sufficiently 
low that a superfluid transition may be observed. In- 
deed, computer simulations yielded evidence of super- 
fluid behavior in very small (less than 20 molecules) P-H2 
clusters (21 and claims have been made of its actual ex- 
perimental observation^ Also, a considerable effort has 



been devoted, in recent times, to the theoretical charac- 
terization of superfluid properties of solvating P-H2 clus- 
ters around linear molecules, such as OCSi^&Si 

The study of hydrogen adsorption on nanocarbons falls 
within the same general research theme, but is also mo- 
tivated by possible practical applications; an important 
example is hydrogen storage, for fueling purposes. So 
far, research along these lines has mostly focused on 
nanotubesj2SiSSi2L22i2£ but it seems worthwhile to extend 
the investigation, possibly providing useful quantitative 
information on adsorption on other nanostructures, in- 
cluding fullerenes. 

In this work, energetic and structural properties of a 
layer of P-H2 molecules adsorbed on a Q fullerene are in- 
vestigated theoretically, by means of ground state Quan- 
tum Monte Carlo (QMC) simulations. In order to pro- 
vide a reasonable, quantitative account of the corrugation 
of the surface of the fullerene, we explicitly modeled in 
our study each individual carbon (C) atom. For compari- 
son, however, we have also performed calculations with a 
simpler model, describing fullerenes as smooth spherical 
surfaces, interacting with P-H2 molecules via an angle- 
averaged potential. 

Only one adsorbed layer is found to be thermodynami- 
cally stable on these small nanocarbons. On a corrugated 
substrate, a commensurate solid layer is observed at equi- 
librium; as the chemical potential is increased, a discon- 
tinuous change to an incommensurate solid layer takes 
place on C20, C36 and Cgo- We could not find, within 
the statistical uncertainties of our calculation, evidence 
of an incommensurate layer on Cso- 

The difference in compression between commensurate 
and incommensurate layers, as measured by the effective 
P-H2 coverage, is approximately 216% for C20, and de- 
creases to ~ 25% for Cgo- 

Obviously, on a smooth fullerene, there is no distinc- 
tion between commensurate and incommensurate layers. 
In the absence of corrugation, energetics of the adsorbed 
layer are determined primarily by the interactions among 
P-H2 molecules. The ground state of p-B-2 in two dimen- 
sions (2D) is a solid, with molecules forming a triangular 
latticeii2£ Our results indicate that p-H 2 molecules at- 
tempt to reproduce the same triangular arrangement as 
on an infinite plane, even when confined to moving on a 
spherical surface of radius as small as a few A. 

Evidence of quantum exchange is absent in these sys- 
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tems, i.e., no evidence suggesting possible supcrfmid be- 
havior is gathered in this work. 

The remainder of this manuscript is organized as fol- 
lows: Sec. |n]offers a description of the model used for our 
system of interest, including a discussion of the potentials 
used and the justifications for underlying assumptions. 
Sec. IIIII involves a brief discussion of the computational 
technique and specific details of its implementation, in 
addition to details of calibration and optimization. The 
results are presented in Sec. II VI finally, Sec. IV! is a 
summary of the findings and our concluding remarks. 



Here, n = 4tt9 a 2 , a being the of the fullercnc and 9 being 
the areal density of C atom on its surface; e and a are 
the parameters of the LJ potential V introduced for the 
fully corrugated model. By using (J2J, one is describing 
the fullerene as a smooth spherical shell, i.e., corrugation 
is neglected. This approximation substantially simpli- 
fies the calculation; it has been adopted in recent studies 
of helium adsorption on buckyballs i ini As mentioned 
above, in this work we have performed calculations based 
on the full model QJ, as well as using the effective poten- 
tial (J2J; results obtained in the two ways are compared 
in Sec. |IV| 



II. MODEL 



We consider a system of N p-H 2 molecules, in the pres- 
ence of a single C; molecule. The latter is assumed fixed 
in space, owing to its relatively large mass; the center of 
the molecule is taken as the origin of a Cartesian coor- 
dinate frame. The I individual C atoms are fixed at po- 
sitions {Rfc}, k=l,2,...,l. All of the atoms and molecules 
are regarded as point particles. The model quantum 
many-body Hamiltonian is therefore the following: 
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Here, m is the mass of a p-H 2 molecule, r« = |r, — r^- 1, 
{rj} (with j=l,2,..,,N) are the positions of the p-H-2 
molecules, V is the potential describing the interaction 
between any two of them, and U represents the interac- 
tion of a P-H2 molecule with a C atom. All pair poten- 
tials are assumed to depend only on relative distances. 
The interaction V is described by the Silvera-Goldman 
potential, which provides an accurate description of en- 
ergetic and structural properties of condensed P-H2 at or- 
dinary conditions of temperature and pressure rfSiSi The 
interaction of a P-H2 molecule and a C atom is modeled 
using a standard 6-12 Lennard-Jones (LJ) potential, with 
e = 32.05 K and a = 3.179 A (see, for instance, Ref. E2). 

The model already contains important physical 
simplifications, such as the neglect of zero-point motion 
of C atoms, as well as the restrictions to additive pairwise 
interactions (to the exclusion of, for example, three-body 
terms), all taken to be central, and the use of the highly 
simplified LJ interaction. On the other hand, is the 
simplest microscopic model that explicitly takes into ac- 
count the corrugation of the surface of the buckyball. 

A further simplification can be introduced by replac- 
ing the third term in Q with JV Ufa), where U is the 
following, spherically symmetric external potential (sec 
Ref. [3 for details): 
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III. COMPUTATIONAL METHOD 

Accurate ground state expectation values for quantum 
many-body systems described by a Hamiltonian such as 
Q can be computed by means of QMC simulations. 
In this work, the method utilized is Variational Path 
Integral (VPI), which is an extension to zero tempera- 
ture of the standard, Path Integral Monte Carlo (PIMC) 
method^ VPI (also referred to as Path Integral Ground 
State, PIGS^i) is a projection technique, which filters 
out the exact ground state wave function out of an initial 
trial state. It is therefore closely related to other ground 
state projection methods, such as Diffusion Monte Carlo 
(DMC), but has a few distinct advantages (for a discus- 
sion, see, for instance, Ref. W^i . 

The VPI technique works as follows: One samples a 
large set {X m }, m = 1,2, ...,M, of discretized many- 
particle paths X = R0R1...R2L, where Rj = rjirj2--.rjN 
represents the positions of all V particles in the system 
at the jth "slice" . These paths are sampled based on the 
probability density 

( 2L-1 , 

V{X)(x^ t (Ro)^t(R2l) J] G(Rj,R J+ i,r) > (3) 

^ j=o J 

where ^t(R) is a trial wave function for the many- 
body system, and G(R,R',t) is a short-time (i.e., 
t — > 0) approximation to the imaginary-time propaga- 
tor (^R cxp — tH R'\. One can showSi that, for any 

choice of V^t, in the limits Lt — > 00, r — > 0, Rl is sam- 
pled from the square of the exact ground state wave func- 
tion $(i?). This allows one to calculate the ground state 
expectation value of any observable O diagonal in the 
coordinate representation as a simple statistical average 
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The ground state energy can be evaluated using the con- 
venient "mixed estimator" , which yields an unbiased es- 
timate: 
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Both expressions Q, (JSJ can be rendered arbitrarily ac- 
curate by letting M — > oo. 

The form of G used here is the primitive approximation 

G{R,R',t)= Pf {R,R',t) e -i^(«)+^(«')] + o( T 3) ) ( 6) 

where £/(i?) is the total potential energy of the system 
at the configuration R, and 



TV 



Pf (R,R',t) = (4tt\t)- 3N ^ [Jexp 
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is the exact propagator for a system of N non-interacting 
particles. 

For a particular choice of t, since the projection time 
A = Lt is necessarily finite, simulations must be per- 
formed for increasing L, until one finds that estimates 
of the observables have converged, within statistical er- 
rors. One must also extrapolate estimates obtained for 
different r, in order to infer results in the r — *■ limit. 
A more accurate form of G (such as that thoroughly dis- 
cussed in Ref. accurate to t 4 ), or a more accurate 
trial wave function ty? allows one to obtain convergence 
with larger r and/or shorter L, but, at least in principle, 
the extrapolated results ought to be independent of the 
choices of G and ^ T m^ In summary, with sufficient com- 
puter time and proper calibration of r and L, one can 
generate estimates of physical observables that are exact 
(errors being only statistical). 

Individual simulations involve generating the set {X m } 
using a random walk through path space. The same path 
sampling techniques of PIMC are adopted; included are 
trial moves involving rigid displacements of an entire par- 
ticle path (displace type moves), and multilevel sampling 
updates (see, for instance, Ref. 13^) . which are divided 
into moves of central path portions (i.e., not including 
an endpoint, namely j=l or j=2L), and moves that up- 
date endpoint configurations (which involve in the 
acceptance test, as do the displace moves). 

The trial wave function utilized is of the Jastrow type: 

, n i 

M=l k=l 

/jj e - u ( r «) 

K i<j 



which fulfills the symmetry requirement of a many- 
boson wave function. In the case of the spherically av- 
eraged model, 1=1 and i?i=0. The pseudo-potentials w 
and u were chosen as follows: 



w(r) 



and u(r) 
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where x = 2 for the spherically averaged potential, and 
x = 5 for the corrugated case. The values of the pa- 
rameters a = 80 A x , 7 = 750 A 5 and x were obtained 
empirically, by minimizing the energy expectation value 
computed in separate variational calculations. 



Using the given *S?t, wc observe convergence of the 
ground state energy estimates with a projection time Lt 
= 0.250 K" 1 , with r = r Q = 1/640 K _1 . Estimates 
for all quantities quoted in this study were obtained us- 
ing t = t , even though convergence of the estimates 
for structural quantities can be typically observed with 
values of r significantly greater than t (for an equal pro- 
jection time). 

VPI calculations for a range of p-H 2 coverages were 
carried out for each of C20, C36, Cqq, and Cgo (we refer 
here to the near spherical isomers of each) . We use an ini- 
tial configuration of para-hydrogen molecules surround- 
ing the fullcrcnc at a distance of approximately twice the 
Q radius (i.e., making sure that no P-H2 start off within 
the buckyball). Because of the strongly attractive char- 
acter of the C;, for a small enough number of surround- 
ing hydrogen molecules the system remains spatially con- 
fined (i.e., P-H2 molecules do not evaporate). Thus, even 
though periodic boundary conditions are used in the sim- 
ulation, they have no effect, so long as a sufficiently large 
simulation cell is used. 

The systematic errors of our calculation are at- 
tributable to the finite projection time Lt and the finite 
time step r. Based on comparisons of results obtained 
from simulations with different values of L and t, we es- 
timate our combined systematic error on the total energy 
per P-H2 molecule to be of the order of 0.6 K or less. 



IV. RESULTS 

Physical quantities of interest include the ground state 
energy per P-H2 molecule e(N) and the radial P-H2 den- 
sity p(r) about the fullerene, as a function of the total 
number N of molecules. Results for e(N) for all fullerenes 
considered in this study are shown in Figs. ^ through 0] 

Let us consider first the case of C20 (shown in Fig. 
[TJ, as it allows us to illustrate some of the general fea- 
tures seen on other fullerenes as well. Filled circles rep- 
resent the energy estimates yielded by the fully corru- 
gated model QJ, whereas filled squares represent esti- 
mates obtained using the spherically averaged potential 
U described in Sec. [n] Solid lines are polynomial fits to 
the numerical data. Also shown in the inset is the num- 
ber of adsorbed molecules TV versus the chemical poten- 
tial fi (in K). This is obtained by first fitting the results 
for e(N), then minimizing the grand canonical potential 
<f>(N) = N(e(N)-fx) with respect to N, for different val- 
ues of p. The chemical potential of bulk solid p-B-2 at zero 
temperature, computed by Quantum Monte Carlo using 
the Silvera-Goldman potential, is /j, = —88 K (from Ref. 
133^ . and that is where all of the N(fi) curves end. 

There are some qualitative similarities, but also impor- 
tant differences (mainly at low N) between the results 
yielded by the two models utilized. Firstly, in the large- 
TV limit they yield essentially the same energy estimates, 
as corrugation becomes unimportant for thick adsorbed 
films. This is the case on all fullerenes. 

Both e(N) curves display a single minimum at a spe- 
cific number N of P-H2 molecules. The minimum cor- 
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FIG. 1: (Color online) Energy per p-Eb molecule e(TV) com- 
puted by VPI, as a function of the number of molecules ad- 
sorbed on a C20 fullerene. Filled squares: results obtained 
with the angularly averaged potential @. Filled circles: re- 
sults obtained by explicitly modeling all carbon atoms in the 
fullerene. Solid lines are polynomial fits to the VPI data. In- 
set shows the number of particles TV plotted as a function of 
the chemical potential /1, for the fully corrugated model. The 
chemical potential of bulk solid p-Eh is —88 K (from Ref. I33T) . 
and that is where the TV(p) curve ends. 
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FIG. 2: (Color online) Same as Fig. Q but for a C36 fullerene. 
The commensurate layer has Nc = 20 p-Eh molecules and is 
compressible, as evidenced by the finite slope of the curve 
TV(/z) for 20 < N < 28. 

responds to the mathematical condition e(TV) = fi(N), 
and physically to the formation of a stable p-H-2 layer. 
On the corrugated C20, a stable layer occurs for a num- 
ber TV = Nc = 12, significantly smaller than that 
(TV = TVs = 22) obtained on a smooth fullerene. More- 
over, the energy per molecule ec{Nc) for the corru- 
gated model, is approximately 10% lower than that on 
a smooth substrate (eg (TVs)), i.e., the corrugated model 
yields stronger P-H2 binding. Qualitatively similar re- 
sults are observed on all fullerenes; however, the differ- 
ence between TVc and TVs is seen to decrease with the 
radius of the fullerene, whereas the difference between 
the minimum energy values in the two models is ~ 10% 
for all fullerenes. 

Corrugation introduces the conceptual distinction be- 
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FIG. 3: (Color online) Same as Fig. but for a Ceo fullerene. 
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FIG. 4: Same as Fig. Q but for a Cgo fullerene. 



tween a layer that is commensurate with the substrate 
(i.e., the outer surface of the fullerene), and one that is 
incommensurate with it. In a commensurate layer, each 
P-H2 molecule sits right on top of the center of one of 
the hexagonal (or, pentagonal) faces of the polyhedron 
formed by C atoms. For example, C20 has the shape of 
a dodecahedron, i.e. it has 12 pentagonal faces, corre- 
sponding to as many adsorption sites for P-H2 molecules. 
In an incommensurate layer, on the other hand, the ar- 
rangement of P-H2 molecules is mostly determined by 
their mutual interactions; the fullerene, in this case, 
merely provides a background attractive potential and 
a curved geometry. 

On all fullerenes studied here, using the corrugated 
model, the equilibrium adsorbed P-H2 layer is found to be 
commensurate. On a corrugated substrate, one normally 
sees a transition from a commensurate to an incommen- 
surate layer, as the chemical potential is increased from 
its equilibrium valued No such transition can be observed 
on a smooth substrate, for which commensuration is un- 
defined. This is indeed what we generally observe on 
the various fullerenes that we have considered. Details, 
however, differ for the different systems. 

On C20, the number of adsorbed P-H2 molecules re- 
mains constant as the chemical potential is increased (i.e., 
no compression of the commensurate layer is observed), 
until it jumps discontinuously from Nc = 12 to TV/ = 22 
(see Fig. QJ. This signals the (abrupt) appearance of 
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0.2021 
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4.11 


42 


0.1979 






45 


0.2120 



TABLE I: Effective p-ID monolayer coverages for both the 
corrugated and uncorrugated models of the fullerenes. For the 
model including corrugation, Nc and Ni mark the number of 
P-H2 molecules in the adsorbed layer at commensuration and 
incommensuration, respectively, with coverages 9c and 9i (in 
A -2 ). The number N$ marks the number of P-H2 molecules 
in the adsorbed layer on a smooth fullerene, 9s being the 
corresponding coverage. The radius a of each fullerene (in A) 
is also given. 
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an incommensurate layer. The number iVj of molecules 
in such an incommensurate layer is the same as in the 
equilibrium layer on a smooth C20; this leads us to sur- 
mise that the physics of the incommensurate layer on the 
corrugated fullerene and that of the equilibrium layer on 
the smooth fullerene are essentially the same. Thus, the 
main effect of the neglect of corrugation associated with 
the spherically averaged potential model, is the absence 
of the commensurate layer. We come back to this point 
later, when discussing density profiles. 

The incommensurate layer that forms on a C20 
fullerene is compressible, i.e., the number of adsorbed 
molecules is seen to increase monotonically with /i, up 
to N ~ 28; no second layer formation is seen before /i 
reaches the value (/i ) corresponding bulk p-H.2, above 
which thicker adsorbed films are thermodynamically un- 
stable. In fact, on none of the fullerenes considered here, 
do we find more than one stable adsorbed layer. 

On C36, our results indicate that the commensurate 
layer is compressible; here too, however, a discontinu- 
ous transition occurs to an incommensurate layer, also 
compressible and again physically similar to the one that 
forms on a smooth fullerene. 

The same physical behavior to that observed on C20 is 
seen on Cgo and Cso, with an incompressible commensu- 
rate layer. On Cso, though, we fail to observe an incom- 
mensurate layer; it should also be noticed that, on this 
system, it is Nc — 42 and N$ = 45, i.e., there is a much 
smaller relative difference between Nc and N$ than on 
the other fullerenes. Our results are summarized in Ta- 
ble U where the numbers are listed of p-tb molecules 
Nc and Ni in the commensurate and incommensurate 
layers adsorbed on a corrugated fullerene, as well as Ns 
for a smooth fullerene. Also shown are the values of the 
effective coverage (2D density) 
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where, again, a is the radius of the fullerene. The effect 
of the curvature of the substrate can be quantitatively 
established by comparing these values of 9 to the equilib- 
rium coverage of P-H2 on a graphite substrate^ estimated 
at 0.070 A -2 . It should be noted that this definition of 
9 is used to compare the coverage as a function of the 



FIG. 5: (Color online) Radial density profiles of a p-ID layer 
adsorbed on Cj. Solid lines: profiles obtained for the smooth 
fullerene model. Dotted lines: profiles of commensurate lay- 
ers on corrugated model. Dashed lines: profiles of incommen- 
surate layers on corrugated models. Density is given in A~ 3 , 
whereas the distance r from the center of the fullerene is given 
in A. 



surface area of the substrate. A definition of the coverage 
incorporating the distance at which the P-H2 sit above 
the fullerenes is given below (the two definitions would 
be equivalent for a planar substrate, a distinction only 
arising with the introduction of non-zero curvature). 

Structural information about these systems is offered 
by the radial density profiles p(r) of P-H2 on each of the 
four Ci for the two models, all shown in Fig. |SJ Solid lines 
denote the radial profiles on a smooth fullerene, while 
dotted and dashed lines denote radial density profiles in 
the corrugated model at commensuration and incommen- 
suration, respectively. These results are qualitatively the 
same for each fullerene. The peak of the profile is shifted 
away from the center of the fullerene, as the system makes 
a transition from a commensurate to an incommensu- 
rate layer, which is not physically unexpected. No ev- 
idence is seen of second layer formation, on any of the 
fullerenes studied here, for values of the chemical poten- 
tial for which the adsorbed film is thermodynamically 
stable (i.e., [i < fi ). 

On comparing the density profile for the incommensu- 
rate layer on a corrugated fullerene (dashed lines) with 
that for a smooth fullerene (solid lines), one observes a 
further shift to the right for the case of a spherically 
averaged potential. Perhaps more interestingly, the in- 
commensurate layer on a corrugated fullerene features a 
greater spatial width, with respect to that on a smooth 
fullerene (particularly on C20 and C36). Otherwise, these 
layers seem physically similar; for example, one can com- 
pute an effective 2D density 9 e ff (different than the cov- 
erage 9 discussed above) defined as 
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where a p is the position of the peak of p(r). The values 
obtained for the two layers are very close. In the case 
of a smooth fullcrenc, as one goes from C20 to Cso, &eff 
approaches from below the value 0.067 A -2 , namely the 
equilibrium density of P-H2 in 2Di2£ This suggests that 
the physics of the incommensurate layer (and that on a 
smooth fullerene) is determined primarily by the interac- 
tion among P-H2 molecules, which attempt to replicate, 
on a curved surface, the same arrangement as in 2D. One 
further thing to note is that the width of the adsorbed 
P-H2 layer on these systems is of the order of 1 A, very 
close to that of an adsorbed monolayer film on graphite. 4 
The computational method adopted here does not al- 
low one to make a direct estimation of the p-ITj molecule 
exchange frequency, unlike its finite temperature counter- 
part (Path Integral Monte Carlo). Nevertheless, from vi- 
sual inspection of many-particle configurations generated 
in the Monte Carlo simulation, we observe little or no 
overlap of paths associated to different molecules, which 
is substantial evidence that many-particle permutations 
are absent in this system. This is consistent with the high 
degree of localization that molecules experience, both in 
the commensurate and in the incommensurate adsorbed 
layers. 

V. CONCLUSIONS 

Using a numerically exact ground state Quantum 
Monte Carlo method, we studied P-H2 adsorption on the 
outer surface of near-spherical fullerenes. We performed 
calculations based on a simple model, in which all car- 
bon atoms are included explicitly, i.e., corrugation of the 
surface of the fullcrenc is captured. 

A single solid layer of P-H2 is found to be thermody- 
namically stable on the fullerenes studied in this work. 
We find that the equilibrium adsorbed layer is commen- 
surate with the corrugated surface of the fullerene. Such 
a layer is found to be compressible on one of the fullerenes 
(namely, C36), incompressible on the others. On increas- 
ing the chemical potential, a discontinuous transition is 
observed to an incommensurate layer on all fullerenes, 



except for the largest one considered here (Cgo). 

The basic physics of the incommensurate layer is driven 
primarily by the interaction among P-H2 molecules, 
which attempt to reproduce the same triangular arrange- 
ment as on an infinite plane, even when confined to mov- 
ing on a spherical surface of radius as small as a few A. 
Indeed, the incommensurate layer is very similar to that 
found using a simpler model of the system, describing 
the fullerene as a smooth spherical substrate. This sim- 
pler model yields results for the energetics and structure 
of the incommensurate layer in good quantitative agree- 
ment with those provided by the fully corrugated model; 
obviously, however, it is necessary to include corrugation 
in order to reproduce the commensurate layer. An in- 
teresting question that arises is whether a commensurate 
layer of helium may exist on these molecules; the theo- 
retical studies performed so far have made use of a spher- 
ically averaged potential to describe the fullerenes^ 

No evidence of quantum exchanges among P-H2 
molecules has been observed in this work. Thus, these 
systems do not appear as not likely candidates for 
the observation of the elusive superfluid phase of para- 
hydrogen. 

In conclusion, this study has provided preliminary in- 
formation on the structure and energetics of hydrogen 
films adsorbed to the exterior of buckyballs, though much 
remains to be answered. The experimental tests of sev- 
eral of the stated predictions could be carried out, for 
example, by measuring the mass of adsorbed hydrogen 
on C/ by examination of excitation spectra in a dipole 
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